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Abstract 

We study the decay widths of the narrow resonances -D*(2317) and Dsj(2460) in the chiral quark 

model, together with the well-known D* and D* mesons. All the parameters in our calculation are 

C^ ' taken from Godfrey and Isgur's quark model except the tt^ — r] mixing angle which is fixed by the 

o' 

CN ' D* decay widths. The calculated electromagnetic decay widths agree with those from other groups 

and the experimental data available quite well. However, the pionic decay widths of -D^j (2317) 

^D I and Dsj{24:60) are too small to fit the experimental data. We suspect that the simple chiral quark 

pion axial- vector interaction Hamiltonian is not suitable for hadron strong decays of Dgj {2317) and 
> 

o 

^ ■ PACS numbers: 12.39.Pn, 13.25.Ft, 13.40.Hq 

O' 

O 

iu: 
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Dsj {2460). 
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I. INTRODUCTION 



The discovery of narrow resonances D*j2317) [l| and Dgj (24:60) [4] raises challenges to the 
quark model. Masses of these two states are about lOOMeV lower than the predictions of the 



a. 



,{*) 



potential quark model [3|. Furthermore, the isospin conserving decay channels D^*' -^ D^*^K 
are forbidden by kinematics. The observed pionic decays into Dg* n^ break the isospin 
symmetry. 

several „o„-co„ve„ti„„al scKe.es. such as .o.ccu.es Q. tet.a,ua.. states UM^, 
or the chiral partners of Dg and D* [9|, |lO| etc., have been proposed (for a review, see 



Refs. 



11 



121 ). But the conventional cs interpretation is still attractive if the experimental 



)(*) 
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141 llq . In the heavy quark 
limit, Dsji23l7) and D,j(2460) naturally form a P-wave doublet J^ = (0+, 1+) with ji = |. 



masses of D).- can be accommodated in the quark model 

7) and 0^^(2460) naturallv form a P-w. 
Couple-channel effects could lead to mass shifts jl6|, |l7[ . The observed pionic decays can 



also be understood thr oug h i] — ir^ mixing from the chiral perturbation theory 
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In a previous work j2i| , we have calculated the pionic decay widths of the D^ ■ mesons 
using the ^Pq stron g de cay model. Another simple decay model to deal pionic decays is the 
chiral quark model [21j. 

In this work, we calculate both the electromagnetic (EM) decays and strong decay widths 
of D*- mesons in the quark model. The meson wave functions are taken from the well-known 
Godfrey and Isgur's model |J| which gives an impressingly good overall description of meson 
states. The chiral couplings of light quarks with tt, rj meson and the isospin violating 







Since the 



ir^-r] mixing parameter are taken from the chiral perturbation theory 
relative momentum is very large in the decays of D*- mesons, the relativistic effect should 
be important. We do not make the non-relativistic reductions of the transition operators. 
Instead we keep Dirac spinors in our calculation. 

In the next section, we present our formalism. The decays of D*- are evaluated and 
compared with experimental data in Sec. III. We give a brief discussion and summary in 
Sec. IV. 



II. DECAYS OF dJ^ in THE QUARK MODEL 

The electromagnetic interaction is standard 

L^ = e^Y^A^ ■ (1) 

In the chiral quark model, the pions interact with the light up and down quarks through 
the axial-vector coupling 

L.. = |f^7M75rvl>-a^0.. (2) 



As pointed out in Refs llSl l22l the reasonable value of g\ ranges from 0.75 to 1. 



The pionic decay of D*- cannot occur from the above isospin conserving interaction since 
Dlj has no u— or d— light flavor quarks in the conventional cs configuration. However, this 
isospin breaking decay can occur through tt^-t] mixing due to the up and down quar k mass 
difference. In the chiral perturbation theory, the rj — -k^ mixing amplitude reads JISL |l9| 

"" ~ 4 m,- '""+'"" ~ ■ ^ ' 

From S'[/i?(3) symmetry, the t] coupling should be 

^. = |f%A^-^'^</'.- (4) 

Thus, we have J3J| 

L^, = —i^e^sYl^sd^T^'' . (5) 

v3/r, 

In our calculation, we do not make the non-relativistic reductions of the transition oper- 
ators. The quark fields are kept in the form with Dirac spinors 

^(^) = E / (0j^ [^^(P' ^)^(P' 'y''^ + kp' 'Mp^ ^)e-^'i ' (6) 

where the anti-commutation relation of quarks reads 

{dip,s),d\p',s')} = {bip,s),b\p',s')} = {2Trf{2p')6s,5{p-p'). (7) 

Accordingly, the meson wave functions are expressed with the quark operators 

1 /" d^k 1 

|M, P) = ^J (^-7^^(^' P^ ^s)?>^QrM4lo) . (8) 



'^^/Pqp: 



Q 





TABLE I: The inner tensor structure matrix of Ds mesons. 




Meson 


JP 


Tm 






1 



Ds 0- 

Dt 1- 

D*j{23l7) 0^ 

Dsji2i60) V 

iPi l^ 

3r l^ 



V2 
a ■ e 

'vf 

cr ■ k 

a ■ ea ■ k 

k ■ e 

ik X cr ■ e 



Here we treat the quark operators in a matrix form 



"i.-f^h'- 4-(!H. (9) 



The inner Qq structure matrix Tm of Dg, D*, D*j{2317) and Dsj(2460) mesons are taken 
from the quark model |2l,|2j|. They are hsted in Table HI The spatial wave functions (p{k, P) 
are normalized as 

In the heavy quark limit, ji = L + Sg is a good quantum number when mg -^ oo. The lowest 
0+ and 1+ excitation states 1)^(2317) and £'^^(2460) have ji = i. They form the (0+, 1+) 
doublet of the P-wave orbital excitation. Dsj{2A60) is an ideal mixture of the ^Pi and ^Pi 
states in this limit 

|D,,(2460),P) = ^\Ds,CPi),P) + sJI\Ds,CP^),P) . (11) 

The spatial wave functions are related to the simple harmonic oscillator (SHO) wave 
functions in Ref. |3| 



^ "fc,P) = 0(fc2). (12) 



For the ground states, 



For the P-wave states. 



3/2 



(fc2) = (h^\ e--''^"' . (13) 



(fc^) = ^(^) e-"W. (14) 



In the rest frame of A, the decay width of a process A ^ B + C is 

T{A^B + C) = -^^ fdQ\^iA^B + C)\\ (15) 

^ B + C) = {BCmo)\A) , (16) 



with 

, I ^ mi - [Ms + McfKMl - {Ms - Mcf)f' 
'^^' 2Ma ■ ^ ^ 

The wave function in Eq. p2j) is calculated in the rest frame, which is an approximation 
valid only for small \P\ for mesons in motion. The calculation of the Lorentz invariant ^ 
matrix element should be done in a suitable frame in which the relativistic effect due to 
|P| is small. In our calculation, the C particle is tt or 7, which is treated as an elementary 
particle. We calculate the invariant matrix element ^ in the frame Pa + Pb = very like 
the Breit frame. If the heavy quark is the spectator, the relevant kinematics are 

^ ' (18) 

(19) 
(20) 
(21) 

(22) 
(23) 
(24) 

(25) 

lltq -r llbQ Z, 

where 

.« = ^^, (26) 

Vq = ^-VQ, (27) 



Pa = 


-,p'c 




Pb = 


-\p'c 




Pq,A = 


k + ^P'c 




Pq,B = 


k - \p'c 




Pq,a = 


-k 




Pq,b = 


-k 




kA = 




, 1 
= k + -VqPc 


ks = 


mQPq,B - rUqPQ^B 


= k- -riQP'c 



and the momentum of C particle in this frame is 



(Ml - {M B + Mcf){Ml-{M B - Mc) 
2M1 + 2Ml - M2 



fC ni\/r2 , r,nf2 _ J\^2 ' y-^"J 

B C 



The matrix elements are listed below. Following Ref. |2J, the formulae are all written 
in the way similar to the non-relativistic formulae except the overlapping integrals which 
approach unity in the non-relativistic limit m —* oo and deviate from unity significantly for 
the light quarks. 

For the pionic decay of the D* meson, we have 

Its radiative decay amphtude reads 

— ^ J = ley -p^xe^ [f^qFsip'^^, m„ tjq) - fiqFsip"^^, tjiq, r]g)] . (30) 

In the case of D*, the following substitution is understood 

D* (2317) can decay into Dg through the emission of one ir. The decay matrix element 
contains two terms 



(P7r|L(0)|5) . [3gAE^(3^. , 



+ ^VQ^^§j-^fF,{pl,m„VQ). (32) 
Its radiative decay matrix elements also contain two pieces 

{V^\L,{0)\S) = (^7^(0)15)^ + {V^\L,{0)\S)m , (33) 

{VjmmS)E _ yi^^^ . ^^ [^^,F,{p'^, m„ vq) - ^,qF,{p% m^, r^,)] , (34) 



\/AEvEs 
{V'i\h-,(a)\S)M flp: 



'?.. 



[ftqVQFaip"^^, ruq, ijq) - fiQijqFsip"^^, mq, r]q)] . (35) 



The decay matrix elements of Dsj{2A60) are complicated, which are listed according to 
its decay modes below. 

• A^V + 71 

^ M/o7T7 F2{p^,mq,r]Q)eyeA 



V^^EvEa V 2 2U m, 



"q 



+ iVQyQ-^^Fi{pl,mq,r]Q)e*yeA. (36) 



A^F + 7 



A 



(P7m(o)iA) 

(P7m(0)|A)^ 
yAEpEA 

(P7m(o)i^)Af 

\/AEpEa 

V^ + 7 



(P7m(o)iA)s+(P7m(o)iA)M, 

W -/?e^ ■ e^ [/ig-i^4(p5, mg, r/g) - ^qF^{p^^, nig, r]g)] 



(37) 
(38) 

• (39) 



(y7m(o)iA) = (F7m(o)iA)i^+ (i/7m(o)iA)M, 

2 






V^EvEa 



'Ip 



2/iQ 



- iVq\l :^-^Fs{p^, niQ, r]g)e*y ■ { p^p^ - -p^l ] ■ ea x e 



1 



A-^S + ^ 



{S^\L,iO)\A) = (^711,(0)1^)^ + {Sj\L,iO)\A)M , 



^/AEaEs 



-(1 + riQ)figF4{p'^^, mg, r]Q) - fiqF^{p% nig, r]Q) 



{Sj\L,{0)\A)m 
v^FJaEs 



+ ^r]qiiQFi{pt^, mg, rig) - HqF^ip"^^, tuq, rig) 

^7 r^2 , , 77 /_2 _ ^ \ I _2 , , 77 /^2 



P'y X S^ * Gy!^ 



(40) 

(41) 



(42) 



(43) 



(44) 



6p' 



■[VQ^qFsip^, rUg, t]q) + rigiigEsip^, tuq, rig)]p^ x e^ ■ Ea ■ (45) 



III. NUMERICAL RESULTS 



In our calculation, the wave function parameter /3 and the quark masses are taken from 
Ref.y 

(3 = 400MeV , rriu = ma = 220MeV , m, = 419MeV , rUc = 1628MeV . (46) 



The value of g\ = 0.87 is taken from Ref. 



2J. We present the pionic decay widths of D'-*^ 



mesons in Table |n] Our results agree with the experimental data very well. 

In Table IIIIl we collect the numerical results of the isospin violating pionic decays of D* 



and D*, states. We used the commonly accepted value 9m = 0.010 for ir'^-r] mixing. The 



TABLE II: Decay widths of D* ^ D + vr in unit of MeV and the relevant Fi values. 



V ^P + t: 


Pn 


Exp. [25] 


i^i 


Present work 


Ref. 24 


D*^ ^ D^+TT^ 


38 


0.030 


0.65 


0.028 


0.029 


D*^ ^ DO + 7r± 


40 


0.065 


0.65 


0.061 


0.064 


D*o ^D^ + TT^ 


43 


< 1.3 


0.65 


0.039 


0.041 



TABLE III: Pionic decay widths of D* and D),- in unit of keV and relevant Fi values. 





P. (MeV) 


Fi 


F2 


Present work 


Dt ^ D, + vrO 


49 


0.80 




7.4 X 10-3 


D*.i2317)^Ds + TT^ 


298 


0.70 


0.52 


1.9 


Dsji2i60) ^D* + ^° 


297 


0.70 


0.52 


1.9 



radiative decay widths are listed in Table HVl together with some results from other groups. 
The relevant overlapping integrals are collected in Table |3 where 

F^' = Fiip'^,m,,VQ) (47) 

F? = Fi{p%mQ,r^,) (48) 

From the table, we see clearly that the overlapping integrals, which are related to the 
relativistic effects, are very important for the light quarks. For the heavy quarks, the 
overlapping integrals always approach unity. 

TABLE IV: Radiative decay widths in unit of keV. 







Present work 


QM [26] QM [27] QM [28] QSR [29] LCSR [30] 


VMD [11] 


£)*± _, D±^ 




0.25 


0.36 


0.050 




0.23 




Z?*o ^ DO7 




14.5 


17.9 


7.3 




12.9 




Dl - D,^ 




0.065 


0.118 


0.101 




0.13 




Di.{2m) ^ 


■Dl-i 


1.5 






1.9 


4-6 


0.85 


i?,j(2460) ^ 


■Dsj 


6.3 






6.2 


19-29 


3.3 


D,j(2460) ^ 


■Dli 


3.7 






5.5 


0.6-1.1 


1.5 


D,j(2460) ^ 


D^(2317) 


0.18 






0.012 


0.5-0.8 





TABLE V: Overlapping integrals related to radiative decays. 







i^l 




F! 


pQ 


F! 


pQ 


Z?*± -^ D^j 




0.39 


0.94 










D*^ -^ D^-f 




0.39 


0.94 










Dl - I?.7 




0.62 


0.94 










D*^.(2317) ^ 


Dll 


0.60 


0.94 


0.56 


0.93 






Z?si(2460) ^ 


Dsl 


0.47 


0.92 


0.44 


0.91 






D,j(2460) ^ 


Dll 


0.54 


0.93 


0.51 


0.92 






i?,j(2460) ^ 


^:, (2317)7 


0.62 


0.94 


0.58 


0.93 


0.50 


0.91 



IV. SUMMARY 

At present, only the decay widths of D* mesons have been measured. For 

D*, Dsj{2317), Dsj{2A60) states, the ratios between their radiative and pionic decay widths 

are available experimentally, which are collected in Table IVII The experimental data are 

taken from "Review of Particle Physics" by Particle Data Group (PDG) 25| and its online 

update server: http://pdg.lbl.gov/pdg.htnil. 

Our calculated ratios — -r ^„ „ and — | ^^ — agree with the experimental ratio 

L)*o -^ DH^ D* -^ Ds-i 

within a factor of two. Such an agreement is quite interesting if we keep in mind both our 
model wave function and the strong decay Hamiltonian are so simple. 

Our calculated ratio -^ -— — ^r is nearly six times smaller than the experimental data. 

This may be partly attributed to the uncertainty of our naive SHO wave functions. In the 
y ^ P + 7 formula Eq. (jHUI) . there exists a strong cancellation between the light and heavy 
quark contributions. We have 

^ ^ T^ ~ _ 1600 ^ 

and 



F? 



0.4 , (50) 



i.e., /UjFg' ^ —jjcF^ which leads to the strong cancellation. The sensitivity of the overlapping 
integrals to the uncertainty of the wave function is amplified in this case. For example, if 
we change the (3 parameter to /? = 300MeV, we have F^ ^ 0.48. Then the resulting ratio 



TABLE VI: Branching ratios between radiative and pionic decays. 



Exp. Present work 



D*.{23n) ^ DsTt'^ 
Dsji2i60) ^ Ds7 

Dsj{2460) -^ D*7r0 
Dsj{2460) ^ D*-f 

Dsj{2A60) ^D^vrO 

i^,j (2460)^^:^.(2317)7 

L>,j(2460) -^ L»*7r0 



0.62 0.35 

0.062 0.11 

0.052 0.009 

< 0.059 0.79 

0.31 3.3 

<0.16 1.9 

< 0.22 0.094 



will increase to ~ 0.02. 

As can be seen in Table IIV| the radiative decay widths of different channels in this work 



are comparable with those from other groups (see also ref. 1311 1. However, the pionic decay 
widths of Dsj{2317) and Dsj{2A60) from the simple chiral quark model in Table UTTl are ten 
times smaller than those from light-cone QCD sum rules approach p^ and the ^Pq decay 
model J2i|. Hence our calculated ratios between EM and pionic decay widths of Dsj{2317) 
and Dsj{2460) mesons are systematically larger than the experimental data by a factor of 
10. Such a large systematic discrepancy cannot easily be ascribed to either the uncertainty 
of the meson wave function or the uncertainty of the value of the 77 — vr" mixing amplitude 



21 



32]. We tend to conclude that the simple strong decay mechanism based on the pion 
and chiral quark axial vector coupling is not realistic if Dsj{2317) and Dsj{2A60) mesons are 
conventional cs states. 

In summary, we perform a systematic calculation of the decay widths of D*, D*(2317), 
and Dsj{2A60). The EM radiative decay widths agree with the available experimental data 
and other model results reasonably well. But the isospin violating pionic decay widths of 
D*A2317) and Dsj{2460) are too small to fit the experimental data. This disagreement 
cannot easily be resolved by changing the wave functions or the rj — n^ mixing amplitude in 
the chiral quark model. One may wonder whether other possible theoretical schemes such 
as the coupled-channel effects, hybrid meson, molecule state or tetraquark interpretations of 



10 



these two resonances may resolve the above discrepancy. However, there is no clear evidence 
in favor of these exotic schemes from BABAR's most recent extensive measurement J33|. 
Therefore we strongly suspect the chiral quark pion interaction Hamiltonian may be too 
simple to describe strong decays reliably. 
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APPENDIX A: THE OVERLAPPING INTEGRALS 



All -Fj(p^, m, r]) can be expressed as 



where 



Fi{p^, 171,7]) = e 4 



r,p2//32 









-fc2//3^ 



' g° + m Iq'^ + m 



2g0 



2q 



/O 



Gi 



Gi 

G2 

Gs 



2Jf 
3/32 



2(fc-p)2-(fc2 + ^p2) 
{q'^ + m){q'^ + m) 

m{q^ + g'° + 2m) 



m{k ■ p) 



(gO + m) (g'O + m) (gO + m) (g'O + m) {q^ + q'°) 

Am{p ■ kY 



m{q^ + q'^ + 2m) 

(gO + m) (g'O + m) {q^ + m) {q'° + m) (g° + q'^) 

m{q° + q'° + 2m) k^ - {k ■ p)^ 

{qO + m){q'0 + m) 'J3^ 

2fc2 r m(g° + q'^ + 2m) 



Am{p ■ ky 



[q^ + m,){q'^ + m) {q'^ + m,){q'^ + m){q^ + q' 



/o\ 



and q and q' are quark's momenta 



q = k + -p 

q = k p 

^ 2^ 
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